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The processes of merging and splitting dilute-gas Bose-Einstein condensates are studied in the 
nonadiabatic, high-density regime. Rich dynamics are found. Depending on the experimental 
parameters, uniform soliton trains containing more than ten solitons or the formation of a high- 
density bulge as well as quantum (or dispersive) shock waves are observed experimentally within 
merged BECs. Our numerical simulations indicate the formation of many vortex rings. In the case 
of splitting a BEC, the transition from sound-wave formation to dispersive shock- wave formation 
is studied by use of increasingly stronger splitting barriers. These experiments realize prototypical 
dispersive shock situations. 
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Dilute-gas Bose-Einstein condensates (BECs) are a 
powerful environment for the study of linear and non- 
linear dynamics. In this context, the merging and split- 
ting of BECs are two particularly important processes 
that, despite their apparent simplicity, provide access to 
a variety of different regimes. For example, the inter- 
ference fringes observed when two BECs collide in the 
low-density regime reveal the fundamental matter wave 
nature underlying BECs [H, 0]. In contrast, numerical 
studies of fast collisions in the high-density regime reveal 
indications of quantum turbulent behavior [3]. Slower 
merging in the high-density regime has experimentally 
been shown to lead to vortex formation 0, Q , providing 
information about the phase evolution of the merging 
BECs. Theoretical studies demonstrate that adiabatic 
merging of BECs can lead to soliton formation, too [6[. 
A precise understanding of the merging dynamics is also 
essential from a technological point of view. For example, 
splitting and merging processes are fundamental opera- 
tions in atom interferometers [7] . The merging of BECs 
is also an important process in the creation of a "con- 
tinuous BEC" [8] where a condensate was continuously 
replenished by newly condensed atoms. 

In this Letter, we study the physics of merging and 
splitting BECs in the nonadiabatic, high-density regime 
and use these processes to investigate dispersive hydro- 
dynamics, including dispersive shock waves (DSWs, also 
referred to as quantum shock waves). After merging two 
BECs, we can observe many solitons (a soliton train). For 
low enough atom numbers, the soliton train is uniform, as 
predicted for a one-dimensional situation 0, E3] • For 
higher atom numbers, a high-density bulge emerges, and 
our numerical simulations suggest that this bulge con- 
sists of many vortex rings due to a transverse instability 
of the soliton train. Splitting a BEC with a repulsive 
barrier can lead to DSWs, and we observe a transition 



from propagating sound waves to DSWs when a suffi- 
ciently strong barrier is used. Finally, we also find shock 
dynamics in a different setting, namely when a high- 
density region in a BEC is suddenly released and allowed 
to spread into a surrounding background of condensed 
atoms. In this case, the initial high-density peak splits 
into two peaks that propagate away from each other, and 
after a sufficient self-steepening of the wavefront, DSWs 
are formed. These are important prototypical situations 
that have been discussed in previous theoretical studies 
ToL 11 . 12, 13, lij], and our results complement previous 



experiments that considered either very narrow initial 
gaps J ii a BEC. produced by a stoppcd-light technique 

ing 



lal . or blast pulses in rotating [17[ and nonrotat- 
14j cylindrical geometries. In addition, we demon- 



strate intriguing similarities between the hydrodynamics 
of BECs and the propagation of laser light through non- 
linear crystals, where related dispersive shock phenomena 
have also been observed [18j . 

All our experiments begin with ultracold clouds 
of 87 Rb atoms in the \F, mp) = 1) hyperfine 
state. The atoms are magnetically contained in an 
elongated Ioffe-Pritchard type trap with frequencies 
{uo x / (27r), uOy Z / (2tt)} = {7, 402} Hz (the x-axis is oriented 
horizontally). Repulsive and attractive barriers for the 
atoms are created with dipole lasers that are far detuned 
from the Rb D-lines at 780 nm and 795 nm. The dipole 
laser beams are sent horizontally through the center of 
the magnetic trap, along the radial (tightly confining) 
^-direction. In the vertical direction (z-axis), the laser 
waist is much larger than the radial extent of the BECs. 
Dynamics are induced in the BEC by rapidly turning a 
dipole beam on or off. To enlarge the resulting features 
we employ an antitrapped expansion 2 ms long before 
imaging [19[. During this expansion, the aspect ratio of 
a BEC formed without the presence of a dipole barrier 
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changes from 57 for the trapped BEC to an aspect ratio 
of about 3. 

In a first set of experiments, the dynamics of merg- 
ing two BECs in the nonadiabatic, high-density regime 
are studied. For this, a dipole beam with a wavelength 
of A = 660 nm, a power of 3.48 mW and waists of 
w x = 27.3 jam and w z = 32.1 jam is used, creating a 
repulsive barrier with a height of 490 nK for the atoms. 
The total atom number for the experiments shown in 
Figs. QJa-f) is about 1 x 10 6 atoms. For a single BEC 
confined in the magnetic trap, this would imply a chem- 
ical potential of \i = 224 nK. Therefore the presence of 
the dipole beam leads to two clearly separated BECs 
(Fig. [U(a)). The beam is turned on before the atoms 
are evaporatively cooled to form a BEC. After a BEC 
has formed on both sides of the barrier and no surround- 
ing thermal cloud is visible, the dipole beam is rapidly 
turned off within less than 250 ns. We let the dynamics 
evolve in the magnetic trap for a variable evolution pe- 
riod before starting the expansion imaging. Directly after 
turning the dipole barrier off, the condensates smoothly 
expand toward each other (Fig. [2(b)). This behavior can 
be described by the well-known dam-breaking problem 
whereby a sharp density gradient develops into a rar- 
efaction wave (as opposed to a shock wave) when the 
background density is zero (see e.g. [14]). Shortly af- 
ter the BECs have collided at the center of the trap, a 
pronounced bulge of higher atom density forms in the 
collision plane (Fig. [2(c)). Subsequently this density 
bulge spreads out from the center of the trap (Figs. [2(c- 
e)). Very pronounced dark notches are observed to form 
within the high-density bulge as shown in Fig. QJd). As 
the bulge spreads out from the collision plane, more and 
more notches are formed to fill the extent of the density 
bulge with an average spacing of roughly 8 jam to 11 /mi. 
After about 55 ms, the bulge and the notches have spread 
over the entire extent of the condensate (Fig. [2(e)). The 
long lifetime, discrete nature, and large amplitude of the 
notches suggests that they are nonlinear coherent struc- 
tures rather than simple sound waves. Our numerical 
simulations show that a soliton train initially develops 
and a bulge region is formed where the solitons decay 
into a large number of vortex rings, see Figs. [2(h-j) and 
(2o| . Experimentally, vortex rings in BECs have been 
observed in [l6|, [2l[ . They are difficult to detect unam- 
biguously in our experimental images that are integrated 
along the line of sight. Fine fringes appear adjacent to 
the bulge region as can be seen, e.g., in Fig. [2(d). The fine 
fringes, together with the steepness of the wave fronts de- 
limiting the density bulge region, are indicative of DSWs. 
The merging process finally results in an axial breathing- 
mode excitation of the BEC. 

The qualitative features of the evolution are fairly in- 
dependent of most experimental parameters. For exam- 
ple, use of two BECs with an initial total atom number 
of 2.2 x 10 6 atoms and a dipole beam with waists of 



FIG. 1: Left column: experimental antitrapped expansion im- 
ages of a BEC collision at t = (a) ms, (b) 5 ms, (c) 10 ms, 
(d) 20 ms, and (e) 55 ms. Right column: (f) integrated cross 
sections of (d). (g) Typical uniform soliton train observed 
for lower atom numbers (experimental image; parameters see 
text), (h-j) Numerical simulations with the same initial pa- 
rameters as in left column. Zoomed-in views after (h) 5 ms, 
(i) 6.25 ms, and (j) 7.5 ms. The simulations show density 
slices of the BEC by the plane z = 0. No antitrapped expan- 
sion was performed in the simulations, see also [201 ]. 



w x = 8.5 /am and w z = 32.1 jam gives qualitatively the 
same results (see Fig. [2Kb)). However, when the atom 
number is strongly reduced, we observe a transition in 
the merging dynamics, both experimentally and numeri- 
cally, from the generation of a high-density bulge to the 
generation of a uniform soliton train with no bulge. A 
typical image of such a soliton train is shown in Fig. [2(g) 
for 22000 atoms, a dipole beam power of 150 /iW, an 
evolution time of 27 ms and an expansion time of 1 ms 
(see also jioj). A one-dimensional analysis akin to that 
in [ll| reveals that the soliton train can be interpreted 
as the result of the interaction of two degenerate rarefac- 
tion waves resulting from the dam-breaking problem. A 
detailed analysis of this transition is beyond the scope of 
this work. 

In a second set of experiments, we investigated the dy- 
namics of splitting a BEC with a repulsive barrier that 
is suddenly turned on in the center of the BEC. For 
very weak barriers that only slightly modify the BEC 
density, the sudden turn-on leads to the propagation of 
sound waves [H, [23[. Strong barriers, in contrast, lead 
to DSWs. We first create BECs with 2.2 x 10 6 atoms 
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FIG. 2: Dynamics induced by repulsive barrier, (a), (b) 
Turning a barrier off after forming a BEC in the presence 
of the barrier, (c), (d) Turning a barrier on after forming 
a BEC without the presence of the barrier, (e), (f) Pulsing 
the barrier on for 1.5 ms after the formation of a BEC with- 
out the presence of the barrier, (a), (c), (e) Weak barrier, 
laser power 360 /xW. (b), (d), (f) Strong barrier, laser power 
1.99 mW. Evolution times before start of imaging procedure: 
(a),(b) 16 ms, (c), (d) 12 ms, (e), (f) 10 ms. 



in the magnetic trap without the presence of the dipole 
barrier. Then a dipole beam with waists of w x = 8.5 /xm 
and w z = 32.1 fim is rapidly turned on and left on for 
a variable evolution time, after which the antitrapped 
expansion procedure is started. The rapid turn-on of 
the dipole beam produces two density peaks that spread 
out to either side, as shown in Fig. [2(c). The measured 
propagation speed of these peaks in the central region of 
the BECs is plotted in Fig. [3] for various powers of the 
dipole beam [24]. For the lowest powers, the speed is 
in full agreement with the calculated longitudinal speed 
of sound, 3.8 mm/s [25[. At these low powers, the den- 
sity peaks are barely visible in the cloud. For stronger 
dipole beams, the propagation speed increases above the 
speed of sound, and for sufficiently strong beams the 
obtained images change qualitatively due to dispersive 
shock formation with two associated speeds. For laser 
powers above approximately 0.6 mW, fine fringes appear 
in front of the propagating peaks. When the dipole beam 
exceeds a power of roughly 1.2 mW, solitons are formed 
in the region between the two wave fronts, as shown in 
Fig. [2(d). For an oblate, rather than cigar-shaped geom- 
etry, related ring-shaped structures have been interpreted 
in the context of DSWs in [3]. Numerical simulations 
[2q( | suggest that this behavior is qualitatively described 
as follows. Two peaks in the density and outward super- 
fluid velocity are generated by the repulsive dipole beam. 
These peaks break when the "entropy" condition is sat- 
isfied [12], causing the generation of two DSWs on the 
inner and outer edges of each peak. Due to a transverse 
instability, these quasi-one-dimensional DSWs break up 
into many vortex rings, leading to interactions. The in- 
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FIG. 3: Speed of wavefront propagation through a BEC. After 
a BEC is formed, a repulsive dipole laser is suddenly turned 
on. The plot shows the propagation speed of the resulting 
wavefronts vs. the applied laser power. Error bars are taken 
from fits of distances vs. time. The line shows the function 
v(P) = a + bVP fitted to the data. Insets show two repre- 
sentative images obtained with an evolution time of 12.5 ms 
and a laser power respectively below and above the power at 
which the wavefronts start breaking into solitons. 



teraction of DSWs has been studied in [11]. Deriving an 
analytic expression for the DSW speeds is complicated 
by the generation of vortex rings and wave interactions. 
Empirically, the dependence of pulse propagation speed 
on the dipole laser power is described reasonably well by 
a square-root dependence, as shown by the fitted curve 
in Fig. [31 

Our experiments reveal that rapidly switching on a 
weak dipole barrier merely leads to two density peaks 
that spread out (Fig. [2(c)), whereas for strong barri- 
ers, soliton formation is observed in the wake behind the 
wavefronts (Fig. [2(d)). A similar transition also exists 
for the merging of BECs: If two BECs are initially sepa- 
rated by a strong barrier, soliton dynamics are observed 
after turning the barrier off, as described in detail above 
(Fig. [T]and Fig. [2(b)). If, however, the initial dipole bar- 
rier is so weak that it merely produces a small density 
suppression in an initial BEC, a density bulge that does 
not contain solitons (Fig. [2(a)) appears after the turn-off. 
Pulsing a barrier on for a short time combines the effects 
of turning on and turning off a barrier (Figs. [2(e), (f)); 
see also [20]. 

Yet another, different prototypical DSW situation 
comes about when a local high-density peak in a BEC 
is suddenly released and is allowed to spread out into a 
surrounding background of lower density. We can realize 
this experimentally by replacing the 660 nm dipole laser 
with an 830 nm laser, leading to an attractive dipole po- 
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(Fig. 3|) followed by the formation of dispersive shock 
waves, are very similar to those observed in nonlinear 
optics. This demonstrates the generality of our results 
and showcases the usefulness of BECs in the study of 
nonlinear wave dynamics. 

PE acknowledges financial support from NSF under 
Grant No. PHY-0652976. 



FIG. 4: Dynamics induced by suddenly turning off a local, 
attractive dipole potential in the center of a BEC. Dipole 
laser wavelength 830 nm, power 61 /xW. Evolution time be- 
tween turn-off and start of expansion: (a) ms, (b) 0.25 ms, 
(c) 1 ms, (d) 14 ms. (e), (f) Development of dispersive shock 
waves when a laser power of 183 fiW is used. Evolution time 
(e) 1 ms, (f) 8 ms. 



tential. In Figs. [U[a-d) such an attractive dipole beam 
with a waist of w x = 5 //m and w z = 41 //m and a 
power of 61 jaW is sent through the center of the mag- 
netic trap in the radial direction. Evaporative cooling in 
the presence of the combined optical and magnetic po- 
tential leads to a BEC with a localized high-density peak 
in its center (Fig. Ufa)). When the dipole beam is sud- 
denly switched off after the formation of the BEC, the 
central high-density peak spreads out into the surround- 
ing parts of the BEC. Interestingly, the high density peak 
does not remain a single peak but quickly splits into two 
peaks (Fig. [4](b)) that subsequently travel outward in op- 
posite directions (Figs. [4fc, d)). Indeed such splitting 
is expected from theory [12[ and has been observed in 
nonlinear optics [181 ] - Following [ll|, one can show that 
in the quasi-one-dimensional regime two counterpropa- 
gating DSWs interacting with trailing rarefaction waves 
are generated in such a situation. For the parameters in 
Figs.[4](a-d), the measured propagation speed of the two 
propagating peaks is 4.34 mm/s. When the initial dipole 
beam power is reduced, the measured speed decreases 
and approaches the longitudinal speed of sound. When 
a stronger dipole beam with a power of 183 /iW is used, 
again two peaks form and spread out (Fig.|H(e)), but now 
DSWs form, marked by solitons in the inner region and 
strong ripples in the outer regions of the BEC (Fig.[4](f)); 
see also [20|. 

In summary, our experiments realize several prototyp- 
ical situations for dispersive shock wave formation and 
demonstrate the transition from sound wave propagation 
to dispersive shock dynamics as increasingly stronger per- 
turbations are applied. Several aspects of this behavior, 
such as the splitting of an initial single peak into two 
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SUPPORTING MATERIAL: NUMERICAL SIMULATIONS 
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FIG. 5: Azimuthally symmetric numerical simulations of the Gross-Pitaevskii equation for the merging experiment. Figs, (a) 
and (b) were generated using the same experimental parameters used to create Fig. 1 in the main text with N — 10 6 particles. 
Fig. (a) is the integrated density. Fig. (b) represents density slices of the condensate by the plane z = 0. This is a zoomed in 
version of the initial interaction process showing complicated spatial structure along with the generation of a large number of 
vortex rings which appears to be due to a transverse (snake) instability of the soliton train. The snake instability of a single 
dark soliton was studied, e.g., in D. L. Feder et al., Phys. Rev. A 62, 053606 (2000). Fig. (c) shows density slices as in (b) and 
involves the same harmonic trapping parameters as in Figs, (a) and (b) but with fewer particles, N — 18, 600 (approximately 
the same as for Fig. 1(g) in the main text). Quasi-one-dimensional behavior is preserved and a soliton train is generated in this 
case. No ant i-t rapped expansion was performed in either simulation. The vertical, y-axis in the figures is measured in units of 
microns. 
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FIG. 6: Fully three-dimensional (x, numerical simulations of the Gross-Pitaevskii equation with the same parameters as 

those used to generate Fig. 1 in the main text but with N = 10 5 particles instead of 10 6 used in the experiment. Fig. (a) 
is the integrated density whereas the Figs, (b) and (c) represent density slices of the condensate by the plane z = and are 
zoomed in relative to (a). The azimuthal symmetry of the solution is visually maintained over the integration times studied. 
A transverse instability of the soliton train and vortex ring generation are clear. No anti-trapped expansion was performed in 
the simulation. The vertical, ^/-axis in the figures is measured in units of microns. 
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FIG. 7: Azimuthally symmetric numerical simulations of the Gross-Pitaevskii equation for the blast wave experiments. Figs, 
(a) and (b) correspond to the same parameters used in the experiment shown in Figs. 2(c) and 2(d) in the main text respectively 
where a repulsive dipole laser was turned on and held on for 12 ms and then turned off for 2 ms. The repulsive dipole beam 
for Fig. (a) has power P — 360/iW whereas the power in (b) is P — 1.99 mW. In the low power case (a), there is no transverse 
instability leading to vortex ring generation and therefore, the excitations are interpreted as sound waves. The dark dots that 
emerge inside the condensates during the evolution in the high power case of Fig. (b) indicate vortex rings due to a transverse 
instability of dispersive shock waves. No anti-trapped expansion was performed in any simulations. The vertical, y-axis in the 
figures is measured in units of microns. 
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FIG. 8: Azimuthally symmetric numerical simulations of the Gross-Pitaevskii equation for the attractive dipole beam exper- 
iment. Fig. (a) corresponds to the same parameters used to create the experimental images in Figs. 4(a-d) in the main text 
(with dipole laser power P = 61 /xW). Fig. (b) corresponds to the same parameters used in the experiments shown in Figs. 4(e) 
and 4(f) (with dipole laser power P = 183 /jlW). In the lower power case (a), there is no transverse instability leading to vortex 
ring generation and therefore, the excitations are interpreted as sound waves. The dark dots that emerge inside the condensates 
during the evolution in the high power case of Fig. (b) indicate vortex rings due to a transverse instability of dispersive shock 
waves. No anti-trapped expansion was performed in any simulations. The vertical, ?/-axis in the figures is measured in units of 
microns. 



